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\^ • Abstract. We present an algorithm for solving stochastic heat equations, whose key 

Mh| ingredient is a non-uniform time discretization of the driving Brownian motion W. For 

this algorithm we derive an error bound in terms of its number of evaluations of one- 

[~^ [ dimensional components of W. The rate of convergence depends on the spatial dimension 

CN ' of the heat equation and on the decay of the eigenfunctions of the covariance of W. 

According to known lower bounds, our algorithm is optimal, up to a constant, and this 

^^ ' optimality cannot be achieved by uniform time discretizations. 

CLi' 

■3 

1. Introduction 

A common technique for the numerical solution of stochastic evolution equations is an 

^ ! Ito-Galerkin approximation, which turns the corresponding infinite-dimensional system of 

Q I stochastic differential equations (SDKs) into a finite-dimensional one. The latter is then 

^ I discretized in time and approximately solved by, e.g., an Euler scheme. 

Q ■ More generally, every numerical algorithm for an evolution equation eventually has to 



O 



^ ■ discretize the driving stochastic process, which frequently is assumed to be a Brownian 

motion on an infinite-dimensional Hilbert space, in space and time. The vast majority of 
algorithms for stochastic evolution equations as well as for SDKs apply a uniform time 
discretization. This means that a finite number of one-dimensional components of the 
driving process are evaluated (simulated) at time instances i/n with a common step-size 
1/n. 

^ ■ In this paper we present and analyze a non-uniform time discretization for a stochastic 

heat equation 



S 



dX{t) = AX{t) dt + B{X{t)) dW{t), 
^^^ X{Q)=i 

on the Hilbert space H = L2(]0, 1[ ). As a key assumption, the system {hi)i(zfqd of eigen- 
functions of the trace class covariance Q of the Brownian motion W coincides with the 
system of eigenfunctions of the Laplace operator A. A finite number of scalar Brownian 
motions {W, hi) is selected, and each of them is evaluated with step-size l/n^ depending 
on its variance. Based on these data, a properly defined implicit Euler scheme is employed 
to compute an approximation X^ to a finite number of components {X, hj) of the solution 
X. Here N denotes the total number of evaluations of scalar Brownian motions {W, hi) 
used by X^^, up to a constant. 
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Our main result is an upper bound for the error 

-1 \\ 1/2 



era=(i5(^ \\Xit)~XUt)rd?j^ 



of X^ in terms of N. The rate of convergence depends on the spatial dimension d and on 
the decay of the eigenvalues of the covariance Q. Assume, for simplicity, that 

LJll-i = \i\2 ' iT'i 

for some 7 G ]d, oo[ \ {2d}, and put 

*. n 1 (2ci-7) + 
« (7,«) = o 



2 2{d + 2) ■ 
Then 

with a constant ci > that only depends on d, 7, B, and ^, see Theorem |2l 

Actually, this upper bound is best possible, not only for the specific algorithm X^ 
but for any algorithm that uses at most a total of A^ evaluations of the scalar Brownian 
motions {W, hi): there exists a constant C2 > that only depends on d, 7, B, and ^ such 
that 

e(X^) > C2 ■ Ar-"*^^''^) 

for any such algorithm. In general, one cannot achieve the optimal rate a*{'y,d) by any 
sequence of algorithms that use a uniform discretization. See Miiller-Gronbach, Ritter 
(2006). 

In the context of stochastic partial differential equations, implicit (Euler) schemes based 
on uniform time discretizations are studied, e.g., by Gyongy (1999), Kloeden, Shott (2001), 
Hausenblas (2002, 2003), Millet, Morien (2005), Walsh (2005), and Yan (2005). Non- 
uniform time discretizations are studied for the first time by Miiller-Gronbach, Ritter 
(2006). In the latter paper, non- uniform time discretizations are used for the numerical 
solution of heat equations with additive noise, i.e., 5 is a function of the time t but not of 
the current value X{t) of the evolution. In this case the solution X is a Gaussian process 
and conditional expectations become feasible as a computational tool. This is no longer 
true for equations with multiplicative noise, as studied in the present paper. Instead, the 
algorithm introduced in the present paper is a general-purpose algorithm. 

Optimality results, as stated here for the algorithm X^, require lower bounds that are 
valid for all (or at least a broad class) of algorithms. For stochastic evolution equations 
the first such lower bound is due to Davie, Gaines (2001), who consider a particular case 
of (^ in spatial dimension d = 1 with a space-time white noise. See Miiller-Gronbach, 
Ritter (2006) for lower bounds for equations (^ in general, with space-time white noise 
for d = 1 and trace class noise for d > 1. 

Our results show the principal significance of non-uniform time discretizations for the 
numerical solution of stochastic evolution equations. Non-uniform and even sequentially 
computed time-discretizations are studied, too, for finite-dimensional systems of SDEs. 
Here, as a rule, those time-discretizations do not improve the order of convergence, but 
only the asymptotic constants. However, improvements may be substantial on the level 
of constants, see Cambanis, Hu (1996), Hofmann, Miiller-Gronbach, Ritter (2001), and 
Miiller-Gronbach (2002, 2004). 
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We outline the content of the paper. In Section |21 we formulate the assumptions on the 
heat equation (|T)) and briefly discuss existence and uniqueness of a mild solution. Our 
algorithm is defined in Section 01 Error bounds and optimality properties are stated in 
Sections m and resp., and proofs are given in Section IHl 

2. Assumptions 

We study stochastic heat equations (|H) on the Hilbert space H = L2(]0, 1[ ). Here 
^ E H for the initial value, and A denotes the Laplace operator with Dirichlet boundary 
conditions on H. Hence Ahi = —fii ■ hi with 

d 

hi{u) = 2''/^-Y[sm{ieTTUe) 
1=1 

and 

2 I -12 
^li = TX ■ |l|2, 

where \i\2 is the Euclidean norm of i G N'^. 

Moreover, W = (iy(t))jg[o_i] denotes a Brownian motion on H, whose covariance Q : 
H ^ H is a trace class operator. Specifically, we assume that Qhi = Xi ■ hi with 

\i = K\i\2) 

for some non-increasing and regularly varying function 

A : [l,oo[->]0,oo[ 

of index —7, where 

7G [d,oo[\{2d} 

and 

/oo 
A(r) ■ r'^~^ dr < 00. 

Note that the latter property always holds if 7 > d. 

Let (■, ■) and || ■ || denote the inner product and the norm, respectively, in H, and consider 
the Hilbert space Hq = Q^^'^(H), equipped with the inner product {Q^^'^hi,Q^^'^h2) = 
(/ii,/i2) for hi,h2 G H. Furthermore, let C = C2{Hq,H) denote the class of Hilbert- 
Schmidt operators from Hq into H, equipped with the Hilbert-Schmidt norm || ■ \\c. We 
assume that the mapping B : H -^ Cis given by pointwise multiplication and a Nemytskij 
operator, i.e., 

ioT X E H and h E Ho, where 

with g G C^(M) such that 



B{x)h = Tg{x) ■ h 
Tg{x) = gox 

\\g'\\oo < 00. 



Remark 1 . We briefly comment on the existence of a mild solution of equation (^ under 
the above conditions on B. 

Note that Tg : H ^ H, see Appell, Zabrejko (1990, Thm. 3.1). Furthermore, Hq C 
Loo(]0, 1[ ), since sup^gjsjd ||/ii||oo < 00, see Manthey, Zausinger (1999, Lemma 2.2). Hence 
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B[x)h E H for every h G Hq. Moreover, (A^ ■ hi)i^^d is a complete orthonormal system 
in Hq and 



J2\\Bix)Xfh,r<snp\\h,\\l-J2\i-\\T, 



9\^) 



|2 



Consequently, -B(x) G £ for every x E H. Clearly, 

(2) ||T,(x)-T,(|/)||<||^'|U-||x-y|| 
for x,y E H, which yields 

(3) \\B{x)-B{y)\\c<K-\\x-y\\ 

with K = supigjsjd ||/ii||oo ■ J2ieNd -^i ' Ik'Hoo- Thus B : H ^ C is Lipschitz continuous. 
Consider the semigroup {S{t))t>o on H that is generated by A, i.e., 

S{t)hi = exp{-fiit) ■ hi. 

From it follows that there exists a continuous process {X{t))teio,i] with values in H, 
which is adapted to the underlying filtration, such that, for every t G [0, 1], 

X{t) = S{t)^+ f S{t-s)B{X{s))dW{s) 
Jo 

holds a.s. This process is uniquely determined a.s., and it is called the mild solution of 
equation ((H). Furthermore, 

(4) sup E\\X{t)f<cu 

*e[o,i] 

where the constant ci > only depends on d,^,X and g. See Da Prato, Zabczyk (1992, 
Sec. 7.1). 

3. The Algorithm 

We construct and analyze an algorithm that is built from the following ingredients: 

(i) an Ito-Galerkin approximation of the stochastic heat equation, 
(ii) a non-uniform time discretization of the corresponding finite-dimensional Brown- 

ian motion, 
(iii) a drift-implicit Euler scheme. 

Put 

(3i{t) = xy'-{w{t),hi) 

for i E N'^ and t G [0,1]. Then (/5i)igi^d is an independent family of standard one- 
dimensional Brownian motions. Let 

Yjit) = {Xit),hj) 

for t G [0, 1] and j G N*^. The real-valued processes Yj = {Yj(t))te[o,i] satisfy the bi-infinite 
system 

dYjit) = -Hj Yj{t) dt+J2 \'^' • {BiXit))hi, hj) d(3i{t) 
Yj{0) = {^,hj) 
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of stochastic differential equations. For any choice of finite sets X, ^7 C N'^ an Ito-Galerkin 
approximation X = {X(t))te[o,i] to X is given by 

with real-valued processes Yj = (^j(t))te[o,i] that solve the finite-dimensional system 

dVjit) = -^ijYjit) dt + J2 ^y^ ■ {B(X{t))hi, hj) d(3i{t) 

(5) iei 

We apply a drift-implicit Euler scheme to the finite-dimensional system ©. This scheme 
is based on a non-uniform discretization of the corresponding finite-dimensional Brownian 
motion, since /3i will be evaluated with step-size l/ui depending on i G X. Accordingly, 
put 

(6) tg^i = i/ui, i = 0,. . .,ni. 

A good choice of the integers rii G N, together with sets X and J', will be presented in 
Section HJ 

In order to understand the construction of this scheme better we first consider a uniform 
discretization, i.e., ti = t^^i = i/n with a common step-size 1/n for all i G X. In this case 
the drift-implicit Euler scheme is given by 

Yj{te) = Yj{te-i) - /x,- F,-(t,) • 1/n + J^ A^' ■ (B[X{U_^))K hj) ■ (Hte) - A(t^-i)) 

iei 

for j E J', where 

(7) x{t) = j2yjit)-hj 

and 

(8) YjiO) = {^,hj). 

Equivalently, 



+ Hh 
In general we define 



iex 



by 



Moreover, we put 



= To < ■ • • < Tm = 1 
{tq, • • • , Tm} = [j{io,i^ ■ ■ ■ , tni,i}- 



iex 



}Cm = {i el : Tm & {to,i, . . . , tni,i}} 

for m = 0, . . . , M, and we define Sm,i for i El and m = 1, . . . , M by 

Sm,i = max({to,i, • • • , tni,i} n [0, r„[). 
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Finally, we use 

1 



(9) r,(t) = nY 



L l + /ij ■ (t Ar^-t Ar^_i) 

for approximation of the semigroup generated by A. Then the drift-implicit Euler scheme 
is given by d?)), (jHl), and 

(10) Yj{t)- ^^^'^ 



^ji^m-l) 






^0^ ■ WiM - Pi{s^,i)) 



for j e J, if 

Equivalently, 

(11) %it)=Tj{t)-{^,hj) 



^ t \'m—lj 'ml • 



iex 



( Yl {B{X{te-,,i))h, hj) ■ ^.^^ ■ (A(t 

\te.i<rm ^^ ^ ^'^^ 



+ E ^i^' ■ E (^(^(^^-M))^i> ^i> ■ F7f^ ■ (f^^it^^i) - A(t.-M)) 



For illustration we consider an example with X = {1, 2}, rii = 6, and ^2 = 4. Then, for 
instance, /C2 = {2}, /C3 = {1}, and /C4 = {1,2}. Moreover, for t G ]r2,r3] = ]l/4, 1/3] the 
approximation X{t) is based on the increments /5i(l/6), /9i(l/3) — /3i(l/6), and /32(l/4), 
while /32(l/2) - /32(l/4) is not used at all. 

4. Error Analysis 

Henceforth constants that are hidden in notations like^ :< and x may only depend on 
d, ^, A and g. 

In the sequel we consider the particular choice 

(12) M - X' 

J = {jen'':\j\2<J} 



in the definition of the approximations X and X. Then the error of the Ito-Galerkin 
approximation X is bounded as follows; see Section |B1 for the proof. 

Proposition 1. For I, J > 

e( [ \\xit)-x{t)\\'dt] ^i/j'+ V Xi/i2i. 

Moreover, we have the following error bound for the implicit Euler scheme X with an 
arbitrary discretization (jH)) specified by n G N"^; again we refer to Section IHl for the proof. 



^Suppose that F and G are functions on some set A with values in [0, 00]. By definition, F{a) < G{a) 
means F{a) < c ■ G{a) for all a G A with some constant c G ]0, oo[. Furthermore, F{a) x G{a) means 
F{a) di G{a) and G{a) ^ F{a). 
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Theorem 1. For I, J > and n eN-^ 

eU \\x{t) - x{t)\\'' dt\ ^i/j'+ Y, v^^+ Y. v-"- 



i\2<I 



\i\2>I 



Suppose that X may use a total of N evaluations of scalar Brownian motions (3i. Then 
a proper choice of / > and n G N"^ is obtained by minimizing 

D{I,n) = ^ \i/ni+ ^ \i/fii 

\i\2<I \i\2>I 

under the constraint ^^i ^jUi < N. Up to a constant, the corresponding optimization 
problem is solved as follows. 
Recall that, by assumption, 

(13) A(r) = r"^ ■ L(r) 

with a slowly varying function L : [1, oo[ -^ ]0, oo[. Let N E N. We take 

to specify the scalar Brownian motions that are evaluated by the algorithm. For i E N'^ 
with \i\2 < I the Brownian motion j3i is evaluated with step-size l/ui where 



rii = riiN 



;,l/2.^1-^.^^(^l/(d+2)))-l/2 



if 7 G [(i, 2d[ and 



rii = riiN 



\]" ■ N 



if 7 G ]2(i, oo[. For the total number of evaluations we thus obtain X]|i|2</^i ^ ^■ 
Moreover, 

inf{D(/,n) : / > 0, n G N^ with ^ n^ < iV} x D{lN,nN) x e^(iV) 

\i\2<I 



with 



^^-7/2 



e,(iV) = N-'/'^^iW . (L(iVi/('^+2)))'/' 

if 7 G [d, 2d[ and 

e,(iV) = iV-V2 

if 7 G ]2d, oo[. See Miiller-Gronbach, Ritter (2006). 
Finally, we take 

J=JN = e:\N). 
Hereby we have completely specified an algorithm X = X^. 



Theorem 2. The error of the algorithm X^ satisfies 

E 



\\X{t)-X%{t)\\'dt]] ^e,{N). 
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Remark 2. The case of a regularly varying functions A of index —2d is not covered by 
Theorem 12 but may be analyzed in a similar way. Assume, for simplicity, that A(r) = r~^'^. 
Take /jv as above, and define 

ni,N= rA]/'-iV/lniVl 

for ieW with \i\2 < In- Note that Y.\i\2<iM ^^-^ ^ ^- ^^^^ 

inf {£)(/, n) : / > 0, n G N^ with ^ n^ < A^} x D{lN,nN) x A^-^ ■ {InNf. 

\i\2<I 

Furthermore, take Jn = N'^'^ ■ (InA^)"^. Due to Theorem Q the resulting algorithm X^ 
satisfies 

eH \\X{t)-X*^{t)fdt\'\ ^N-^/^-\nN. 

Remark 3. Consider the implicit Euler scheme X with a uniform time discretization ©, 
i.e., ti = te^i = i/n for all i G N'' with \i\2 < I and some constant n G N. Assume, for 
simplicity, that A(r) = r^'^ with 7 G ]d, cxd[. By Theorem Q 

e( f \\X{t) - X{t) f dt\ < 1/ J2 + ci(/, n) 
with 

|i|2</ |i|2>/ 

Minimization of this quantity, up to a constant, subject to the constraint n ■ #T < A^ 
leads to 

mi{d{I, n) : I >0, n G N with n ■ #J < A^} x d(/^, njy) x A^ ^7+2 
with 

and 

Take 

J = JjV = iV'/'"2(^T2)^ 

and let X]^"^' denote the resulting algorithm. By definition, un ■ H^In x A^ for the total 
number of evaluations of scalar Brownian motions (3i used by X^\ and Theorem [T] yields 

(14) (e(J \\X{t) - X'}^\t)\\'' dt\\ ^ Ar"'/'+2(^T2) . 

Remark 4. We compare the implicit Euler schemes X^j and X™', both of which roughly 
use N evaluations of scalar Brownian motions. Assume that A(r) = r~"' with 7 G ](i, oo[ \ 
{2(i}, and put 



2 2(rf + 2) ' '" ' 2 2(7 + 2) 
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From Theorem |21 and Remark El we get 

•1 \ \ 1/2 



eH \\X{t)-XUt)fdt\'\ ^iV-"*(^" 



for the non-uniform discretization, and 

»1 \ \ 1/2 



eH \\X{t)-X^''\t)fdt\\ ^AT-^^T'-^) 



for the uniform discretization. 
We always have 

a;*(7, d) > a{-y, d). 

In the hmit for a low degree of smoothness 

lim a*(7, d) = lim a{'y, d) = l/{d + 2). 

7— >d-|- 7^d+ 

Conversely, for a high degree of smoothness 

lim a(7, d) = 1/2, 

7— ►oo 

while a* (7, d) = 1/2 already holds if 7 > 2d. 

5. Optimality 

The results from Section 0] provide upper bounds for the error of specific algorithms. 
In particular, the comparison of the implicit Euler schemes based on uniform and non- 
uniform time discretizations is in fact a comparison only of the corresponding upper 
bounds. It is therefore important to know whether these upper bounds are lower bounds 
for the error as well, and, even more, to raise the following questions: 

(i) Does there exist any algorithm X^r that uses a total of N evaluations of scalar 
Brownian motions Pi and achieves an error significantly smaller than the upper 
bound e*(A^) for the algorithm X^ff 
(ii) Are non-uniform time discretizations superior to uniform ones? 

To answer these questions we consider arbitrary methods that evaluate a finite number 
of Brownian motions jSi at a finite number of points and then produce a curve in H that 
is close to the corresponding realization of X. In general, the selection and evaluation of 
the scalar Brownian motions /5i, is specified by a finite set 

and nodes 

< ti,i < ■ ■ ■ < t„,,i < 1 

for i E I and rii G N. Every Brownian motion j3i with i G X is evaluated at the corre- 
sponding nodes ti^i. The total number of evaluations is given by 

|n|i = y^n|. 
iei 

Formally, an approximation X to X is given by 

where 

(j) : [0, 1] X mI"'i -^ H 
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is any measurable mapping and X = {ii, . . . , ij^}. Here may depend in any way on the 
initial value ^, the eigenvalues A^, and the function g, which is used to define the mapping 
B in the heat equation (P). The error of X is defined by 

1/2 



e{X)=(^E(^J^ \\Xit)-X{t)rd?j^ 



cf . Theorems ^ and El and the subsequent Remarks. 

Let Xn denote the class of all algorithms (fT3|) that use a total of A^ evaluations of the 
scalar Brownian motions f3i, i.e., \n\i = N. We wish to minimize the error in this class, 
and hence we study the A^th minimal error 

e(A^) = inf e{X). 

XgXn 

In particular, our algorithm X^ is of the form ()15|) . and its total number of evaluations 
of scalar Brownian motions is roughly given by A^. 
We obtain a negative answer to Question (i). 

Theorem 3. The sequence of algorithms X^ is asymptotically optimal, i.e., 

e{X*^) X e{N), 

and 

e{N) X e,{N). 

Proof. In view of Theorem |21 it remains to show that e(A^) >z e*(A^), and this lower bound 
is is a consequence of a more general result established in Miiller-Gronbach, Ritter (2006, 
Thm. 1). D 

With respect to Question (ii) one needs to study the subclass X^' C Xtv of algorithms 
that are based on a uniform discretization, i.e., t^^i = i/n for all i & I and some constant 
n G N. The corresponding A^th minimal error in this class is given by 

e"°'(Ar) = inf e(X). 

Remark 5. Consider the specific equation 

dX{t) = AX{t)dt + dW{t), 
^^'^ X(0) = 0, 

i.e., g = 1 or, equivalently, B{x) = id, and assume that A(r) = r^'' with 7 G ]d, oo[\ {2d}. 
Then 

(17) e""(Ar) y iV^^/^+WTs) 

see Miiller-Gronbach, Ritter (2006, Remark 6), so that 

e(X]^") X e^-^X^) X N'^^^^2i^TT) 

follows from Remark El 

We thus conclude that the upper bound (|TH) for the error of the implicit Euler scheme 
X^'^ is sharp. Moreover, these algorithms form an asymptotically optimal sequence among 
all algorithms that use uniform time discretizations. Our comparison of orders of conver- 
gence in Remark m is therefore a result on minimal errors and clearly shows the superiority 
of non-uniform time discretizations. 
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Note that the conclusions from the previous remark only apply to the specific equation 
()16|) . We conjecture, however, that the lower bound (fTTj) holds in general, in which case 
these conclusions hold in general as well. 

6. Proofs 

We start with regularity properties of the solution X of equation (^. For the mean- 
square smoothness of X we have 

(18) E||X(s)-X(t)f ^ |t-s| ■(l + ?/'(min(s,t))) 

for s, t e [0, 1] with 

7A(t) = 5^/i,-i?((X(t),/i,)2) 

satisfying 

(19) / V'(t) dt < oo. 

Jo 

See Miiller-Gronbach, Ritter (2006, Lemma 1). 

Consider the Sobolev space W2 = W2QO, 1[ ) and its subspace W2' ■ Note that hi G 

W2'° and 

(x, hi)ivi = {1 + fii) ■ {x, hi) 

for every x E W2' ■ Consequently, the functions (l + /ii)~^/^ -hi form a complete orthonor- 
mal system in W2' , and 

(20) 1^2 '° = {xeH -.J^l^i-i^^ hi? < 00} 

as well as 

(21) \\x\\l.,<2■Y,^^i■{x^hi)^ 



2 



2 



for X G W2 ' , which is the Poincare inequality. 
Lemma 1. For Lebesgue-almost every t G [0, 1] we have 

X{t) G 1^2'° 
with probability one and 

J2 ^/N ■ E{Tg{X{t)) ■ hi, hjf < l//ii ■ (1 + E\\X 

for every i eN'^. Moreover, 

1 

E\\X{t)\\l,^dt< 00. 
^ 

Proof. Combine (jl9|) . (|20p . and (J2H) to obtain the first and the last claim. 
For the proof of the second claim we note that 

(22) \\Tg{x)\\w^^l + \\x\\w^ 

for X G W2, see Appell, Zabrejko (1990, Theorems 9.2 and 9.5). Furthermore, we may 
assume g{0) = without loss of generality. Then 

(23) T,{W^^') C 1^2''° 
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is easily verified. In view of ()22|) . ()23|) and tlie first statement in the lemma it suffices to 
show that 



(24) 5^1//i,--(x-/i,,/i,f ^l//i, 

for all X e Wl'^ and i G W. 



•^ T/t/1 



To this end fix i, j G N and ^ G {1, . . . , d}, and put 



/i = l/{ten) ■ £-hi, fj = l/ijiTc) ■ £-hj. 



Then 



h ■ \x, hi ■ hj) ^ [x, jr-ji ■ hj) 



dueJ-^ '"31 



^{£xji-hj)'+jl{xji-fj)\ 



Hereby 



zj, 



+ \\X- f:" 



t\\ 1 



and we conclude that 

Hi ■ ^ l//ij • {x ■ hi, hj)^ ^ "™"^ 

which yields i^^. D 



Flliyi' 



6.1. Properties of the Ito-Galerkin approximation. Let Pj and Pj denote the or- 
thogonal projections onto the subspaces spanj/i^ : i G X} and span{hj : j G J'}, respec- 
tively, and put 

B{x) = Pj o B{x) o Pj. 

Then B : H -^ C satisfies (jH)) and X is the mild solution of (^ with B being replaced by 
B. Hence 

(25) sup E\\X{t)f < ci, 

<e[o,i] 

see dH). _ 

We establish an error bound for piecewise constant interpolation of X. 

Lemma 2. For l^JdW and m ^^ 

J2 E\\X{t)-X{i/m)fdt^l/m. 

Proof. Note that ()18j) and ()19|1 are valid, too, for X and 

V^(t) = ^/i,--E(F,'(t)) 

3(^J 

instead of X and ■?/', respectively. For £ E {0, . . . , m — 1} take si G [£/m, (£ + l)/m] with 



i){si)/m < / ?/'(t) dt. 

Je/m 
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On the first subinterval, 

/ E\\X{t)-X{0)fdt<2/m- sup E||X(t)f ^ 1/m, 
Jo te[o,i] 

see (j25|) . On the subintervals [i/m, {i + l)/m] with £ > 1 we proceed as follows. If t G 
[i/m, Si], then 

E\\X{t) -X{i/m)\\' ^ E\\X{t) -X{se-iW + E\\X{se-i) -X{i/m)f 

:< 1/m + / il){s) ds. 

J{e-l)/m 

lit e [se,{i+ l)/m], then 
E\\X{t) -X{i/m)r 

^ E\\X{t) - X{s,)f + E\\X{s,) - X(s,_i)f + E\\X{se-i) - X{£/m)f 

M£+l)/m_ 

:< 1/m + / 'ijj{s) ds. 

J(e-l)/m 

We conclude that 

/•(£+l)/m i-{i+l)/m_ 

/ E\\X{t)-X{i/m)\\'^dt^l/m^ + l/m- tlj{s)ds, 

which completes the proof. D 

Proof of Proposition Recall the particular choice (J12j) of the sets X and JT and let c 
denote the right-hand side in Proposition ^ Moreover, let 

|j|2<J 

for k = 1,2 with 

^■'^ W = E ^y' ■ [ ^M-Hit - ^)) ■ {UX{s)) ■ h,, hj) d(3i{s) 

\i\2>I •'^ 

and 

Y^'\t) = exp(-/i,-t) ■ {i, hj) + J2 >^f ■ [ exp(-/i,-(t - s)) ■ (T,(X(s)) • h,, hj) d^i{s). 

Then 

x(t) = x«(t) + x(2)(t)+ Y, Ut)-hj, 

\3\2>J 

and consequently 

E\\X{s) -X{s)f ds 

- J2 [ E{Y}{t))dt+ [ E\\X^^\t)f dt + [ E\\X^^\s)-X{s)fds. 

\j\2>J^^ Jo Jo 
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We have 

/ E{YJ'\t)Ydt:< J2 ^ilH- I E{Tg{X{t))-hi,hjfdt, 

and therefore 

f' E\\x^'\t)r dt < J2 f E{Y^\t)fdt^ Y. v/ii<c 

by Lemma H Furthermore, 

E f E{Yf{t))dt<llJ'<c 

follows from (jSj), (jlj), and sup^^j^d ||/ii||oo < oo. Finally, if \j\2 < J, then 

E{YJ'\t)-Yj{t)Y< 5^ A,- f E{T,{X{s)) -T,(X{s)),h,hj)' ds, 

m2<-i 

and due to ^ 

E\\X^'\t) -X{t)f ^ f E\\Tg{X{s))-Tg{X{s))fds ^ f E\\X{s) -X{s)f ds 

Jo Jo 

^2c+ [ E\\X^^\s) -X{s)f ds. 
Jo 

Since E(X«(t),X(2)(t)) = 0, we get E\\X^^^\t)f < E\\X{t)f. Use OH) and ^ to 
conclude that 

sup E||X(2)(t)-X(t)f <oo. 
te[o,i] 

It remains to apply Gronwall's Lemma to complete the proof. D 

6.2. Properties of the implicit Euler scheme. Recall the definition © of F^ used 
for approximation of the semigroup. 

Lemma 3. Suppose that i E I and j G J^. Then, for i = 0, . . . ,ni — 1, 

^' F^(t) 



r^ L ■it) 

/ -^dt<2/f,j 



as well as 

■^ f r,-(t) 



/ f^^4rT-exp(-/ii(t-t,,i))') dt^l/n* 
J hi \^j\tt,i) J 



where 

n* = maxjrii : * G X}. 

Furthermore, for < s < t < 1, 



1 r,-(t) 



Tjis) 



< min(l,/ij ■ (t — s)). 
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Proof. For t G [tk,i, tfc+i.J with k > i 

r.-(t) ^ 1 



and therefore 



1 Y^.{t) ^ "^-' 



Thus, if iJij/rii > 1 



1 r2 



3 

and otherwise 









dt<---. TTTT^^rr? < 2//ii, 



too. 

For the proof of the second statement put k* = [t^ ,j ■ n*] and 

T-(t) 
f{t) = -j)^-exp{-fxj{t-te,i)). 

Then < / < 1 and 

>1 "*-l /•(k+l)/n 






It remains to show that 

n*-l 

(26) V sup f{t)^l. 

k^k* te[k/n*,{k+l)/n*] 

To this end assume that t G [k/n*, {k + !)/'«.*] for some k > k* in the sequeL Use 

r,-(t) ^1 1 1 



Tj{te,i) - l + l^j- {k*/n* - t^,^) (1 + fij/n*)''-''* 1 + /x,- ■ (t - k/n*) 
to obtain 

fit) < ((TT^^^y^I^jfeZF - exp(-/ii(A; - r)K)^ + eM-Nik - r)/n*) ■ (/o + /i(t)) 

with 

/o = — 7777-^ — -7 - exp{-fij{k*/n* - te,i)) 

1 + /ij ■ {k*/n* - t^,i) 

and 

^^^'^ = l + ^j.{t-k/n*) - '"P^-^^-^' - ^/"*))- 
Note that 1/(1 + -u) — exp(— -u) ^ min(l/(l + u),u'^) for -u > 0. Let -u = fij/n*. Then 

7- ; — , , ,, — exT){—aAk — k*)/n*) ^ -:; -^ — -—- ■ min(l/(l + u),u^) 

{1 + lij/n*)k-k ^^ ^^^ " ' - {\^u)^~^-^ ^ '^ " ' 

for k > k*, and therefore 



n*-l . ^ 

^ 1(1 + «)'=-'=* 



exp(-(A; -k*).u)) ^ .\^^J ^ ■ min(l/(l + u)\ u') ^ 1. 
k=k' — -- / (l + «)2-l 
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Since max(/o, /i(t)) :< min(l,'u^), we have 



n* — 1 ^ 

J2 exp(-2/x,-(A; - k*)/n*) ■ (/o + h{t)) ^ — 

k=k* ^^ 



which completes the proof of ()26|) . 

For the proof of the third statement let s < t, and assume that s G [T^-ijTf^] and 
t E [Ty-i,Ty]. By definition 

s - r^-i) T-T 1 1 



TjW _ l + /ij- (s-r«-i 



which implies 

1 +/ii ■ (tk -r^-i^ 



r,.(.i 



I/-1 



+ i: 



(, = K+1 



1 



1 



1 + yUj ■ (t-^ - r,_ 



1 



1 + lij- {t- r^_ 



</ij-(t-s). 
Finally, < Tj{t)/Tj{s) < 1. 

Put 

aij{t)=E[(Tg{X{t))-h,,hjY 

and note that 

(27) 2]]ai,,-(t) ^1 + E||X^^^"2 

due to © and sup^^j^d ||/ii||oo ^ 1- 
Lemma 4. For J, J' C N'^ and n G N-^ 



D 



sup Fi|X(t)r ^1. 
te[o,i] 

Proof. At first we slightly modify the process X by replacing the Brownian increments 
Piijm) — Pi{sm,i) in the definition (fTUI) of K,- by increments (3i(t) — Pi{sm,i)- More precisely, 
we consider X{t) = J2jej '^ji^) ' ^j with Yj{0) = (^, hj) and 



Yj{t) 



^ji'Tm-l) 



\ ieJCm 






iPiit) - PiiSm,i)) 



for t G ]Tm_i,Tm,]. Note that Yj and 1^- as well as X and X coincide at the points r„ 
Moreover, by construction of these processes we have 

(28) Yj{Tm) and X^Tm) are measurable w.r.t. a {{Pi{te,i) : t^,i < r^, * G X}) . 
We claim that 

(29) sup E||X(t)f ^1. 

te[o,i] 
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Assume that t G ]Tm-i,Tm] in the following. Observing (|^ we obtain 






< E{Y?{rm-i)) + E Ai/rii ■ aij{sm,i)- 



From (j27p we therefore get 

E||X(t) -X(r„_i)f ^ 1 + max ^||X(r,)|p 

fc=0,...,m— 1 

and we conclude that 

f{s)= sup i?||X(r)||2 

re[0,s] 

is finite for s G [0, 1], since E\\X{0)f = \\^f < oo. 
Analogously to (fTTj) we have 

Yjit) = Tjit).{^,hj) 
which implies 






ieX \tl.i<Tm J '* 

due to the measurability property (J2H1)- Use (j77j) to derive 

i?iix(t)f ^iieir+E^i- ( E (n-/(^^~M))-(^At,,,-t,_M) 

Jo 
so that (j29|) follows by means of Gronwall's Lemma. 

For the process X we apply ()11|) and observe ()28p again to obtain 

E aiAti-i,i) ■ fT^. 
Using (J29|) we conclude that 

We turn to the mean-square regularity of the process X. 



D 
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Lemma 5. For X,JcN'^,nEN^, andO<s<t<l 

E\\X{s) - X{t)f ^ (t - s) ■ (1 + i^is)) + Y, ^i/^i^ 

iei 

where 

3&J 

Moreover, 

(31) / 4i{s)ds^ 1. 



^0 

Proof. Since s G Jt^-i, t^] and tii < r^ implies t^_i,i < s, we obtain 



3 

rii-l 



Jo Jo iei V£=o -/t£,i-^7l^ 

(rii-l 
ieX £=0 

from ()30p and Lemma El It follows that 

rii — l 

i-Mds d Kf- + E V". ■ i](i + B(l|.?fe)f )) 



' ^^<-'' ., 



iei i=o 

see (|77|l . Use Lemma El to complete the proof of (jHT|) . 

Assume that s < t with s G [Tm_i, Tm] and t G ]tk_i, r^] for m < k. Then 

where 

/Ci(s, t) = {£ G {1, . . . , rii} : ti^i G Jr^, r«;]} 

if s > Tjn-i and 

/Ci(s, t) = {^ G {1, . . . , Hi} : t£,i G [r™, r^]} 

if s = Tm-i- By Lemma El 

EiYjis) - F,-(t))2 ^ /x,- -it-s)- EiY^{s)) + J] A^/n, ■ I ^ «^,.-(^^-i,i) 

iex \i(^iCi(s,t) 

Note that #/Ci(s,t) < 1 + n^ ■ (t — s), and apply (P7|) together with Lemma JH to obtain 

E\\X{s) - X(t)f <{t-s)- i;{s) + X ^iK ■ #^i(^. ^) 






as claimed. D 
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In view of Lemma |3 and Lemma El we may proceed as in the proof of Lemma ^ to 
obtain the following error bound for piecewise constant interpolation of X. 

Lemma 6. For l,J<zW,ne N^, and i el 

J2 E\\X{t)-X{te-iA)fdt^l/ni + J2>^i'/^i'- 

Proof of Theorem^ Recall the particular choice p2|) of the sets X and J and consider 
the corresponding Ito-Galerkin approximation X. Because of Proposition [T] it suffices to 
show that 

(32) j E\\X{t)-X{t)\\'^dt^ Y^ Xi/m. 

For z/ = 1, 2, 3 we define 



Ut\t) = E ^\'' ■ I E ^fi/(^'^) ■ !]*..*..,.] (^) df^^(') 



with 



\i\2<I ° i=0 



V;« (., t) = exp(-/.,- ■ {t - s)) ■ {Tg{X{s)) - T,(X(t,,,)) • hi, hj), 
vS/s, t) = eM-H ■ it - s)) ■ {T,(Xite,i)) - r,(X(t,,,)) • h^, hj), 
K3) ,„ .^ _ /„.., „ ,. „^^ ^jit) 



V,^;//s,t)=\^exp{-f^j.{t-s))~.^^j-{T,{X{te,i))-h,hj). 
Furthermore, we put 

Ufit)= E Af ■fT^-^^^^^^^'".^))-^-^^-)-^^^^^)-^^^^-'^)) 

and 

Ufit) = E A]/^ ■ :p^/^ ■ (T,(X(s„,,)) ■ h, hj) . {f3i{rj - A(t)) 

if t G ]Tm-i, Tm]. Then, by definition, 

= (eM-Ht) - ri(t)) ■ (e, h^) + f/f (t) + uf{t) + uf\t) + uf{t) - uf\t). 

We separately estimate the terms from the right-hand side of this equation. 
Lemma El yields 

(33) E (<^' ^i)' • / i^M-Ht) - ^At)? dt) ^ 1/n* ^ E ^^Z"- 

\3\2<J V >^0 / |.|^<^ 

By (0) and Lemma H 

(34) E EiU^'\t))"^ E ^^- E / E\\X{s)-X{t,,i)rds]^ E Vni. 

Put 

f{s)=E\\X{s)-X{s)r, 
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which is finite because of ()25p and Lemma |3] By Q, Lemma |2l and Lemma IHl 

(35) 

lil2<J 

|i|2</ V £=0 ■^*^*«.'^ 

^ E v^^+ / ^(^)^^- 

|i|2</ "^° 

Suppose that s G ]t^,i, t^+i,i]. Then 

|exp(-/ij(t - s)) - exp(-/ij(t - t£_i))| < exp(-/ij(t - s)) ■ /ij/n^ 
and therefore 

/^(exp(-,,.(t-.))-^)^t 

follows from Lemma El Hereby 

E{uf\t)fdt 



^ , , 1 /-t^i-l / -p /+A \ 2 

-< 






K|2</ \ ° ^=0 

which implies 

(36) E / E{uf{t)ydt^ E ^^/"^ 

see (j27|) and Lemma |3] 
By the same facts, 

(37) E E{uf(t))"^ E V^i 

\3\2<J \i\2<I 

and 

(38) E E{uf{t)r ^ E ^^/^- 

|j|2<J |i|2</ 

if t G ]T:;„-i,rm]. 
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Combining (|HH j) - (jnHj) we obtain 

pr pr pt 

/ f{t)dt< ^\i/ni+ / / f{s)dsdt. 
Jo |.|^<^ Jo Jo 

Finally, apply Gronwall's Lemma to derive Jg f{t) dt ^ X]|i|2</ ^i/^i^ &s claimed in (jH2I)- 

D 
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